Switching energy-delay of all spin logic devices 
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The need to find low power alternatives to digital electronic circuits has led to increasing interest 
in alternative switching schemes like the magnetic quantum cellular automata(MQCA) that store 
information in nanomagnets which communicate through their magnetic fields. A recent proposal 
called all spin logic (ASL) proposes to communicate between nanomagnets using spin currents which 
are spatially localized and can be conveniently routed. The objective of this paper is to present a 
model for ASL devices that is based on established physics and is benchmarked against available 
experimental data and to use it to investigate switching energy-delay of ASL devices. 



Digital electronic circuits store information in the 
form of capacitor charges that are manipulated using 
transistor-based switches. Switches of this type cur- 
rently operate with a supply voltage of one volt involv- 
ing ^ 10^ — 10^ electrons, requiring 1 — 10 femto- Joules 
(fJs), dissipating 1-10 jiW per switch if operating at 1 
GHz. This dissipation per switch is believed to be the 
single most important impediment to continued minia- 
turization and there is a serious attempt to "reinvent the 
transistor"^ so as to operate at lower voltages. 

A more radical approach is to replace the entire charge- 
based architecture with an architecture based on some 
other state variable such as spin^. For example, MQCA^ 
uses nanomagnets to represent digital information (0 
and 1). Recently an all spin logic (ASL) device^ has 
been proposed whereby information is similarly stored in 
nanomagnets but is communicated via spin currents that 
are spatially localized and can be conveniently routed 
within a spin- coherence length which can be lOO's of 
nanometers^ to micronJ^. 

It has been argued that ASL devices could potentially 
lead to ultralow power switches since a stable nanomag- 
net with an activation barrier of 40 kT could be switched 
with less than an atto Joule (aJ)^. Experimentally, how- 
ever, nanomagnet memory devices typically require tens 
of fJs to switch at speeds that are a factor of 100 to 1000 
lower, raising questions about the potential of ASL de- 
vices to provide a low-power alternative to today's tran- 
sistors. This is because most of the dissipation in switch- 
ing magnets is associated not with the dynamics of mag- 
nets but with the spin transport process and we need a 
suitable model that incorporates both to make reliable 
predictions. This paper presents such a model that is 
based on established physics and is benchmarked against 
the recent experimental result of Yang et al.^. 

In general, the switching energy and energy-delay can 
be written as: 
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V and / are the charge voltage and current respectively 
and tsw is the switching delay. Qtot = I^sw is the total 
charge involved in a switching event. Equation 1 per- 
mits a simple comparison with charge-based devices like 




FIG. 1: (a) An ASL device consisted of input and output 
magnets, (b) Illustrates the self-consistent model, (c) Shows 
the conductance matrices describing the spin-transport. 



today's transistors where Qtot is the amount of charge be- 
ing switched. ASL devices permit low voltage operation 
especially if metallic channels are used. For example, the 
experiment in Ref.[7] uses a switching voltage of ~ 30 
mV and it is of the same order or less for GMR devices, 
far lower than today's transistors, and one objective of 
this paper is to use our quantitative model to provide 
insight into the factors that determine Qtot- 

A generic ASL device is shown in Fig{l^, with charge 
current going through an input magnet (mi) andean ac- 
companying spin-current resulting in spin-torque^ which 
if large enough could flip the output magnet (m2). An- 
alyzing such a device involves coupling (Figjl]3) a model 
for magnetization dynamics described using the Landau- 
Lifshitz-Gilbert (LLG) equation with a spin transport 
model (Figjl]^). For the latter, we adopt what we 
could call a "spin-circuit" approach by combining (a) 
the well-established spin-diffusion model developed by 
Johnson-Silsbe^ and Valet-FertP^l that are now widely 
usecP for spin transport in long channels, with (b) a 
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conductance model for the channel-magnet interface pi- 
oneered by Brataas et. al^, whereby a 4-component 
voltage drop is related to a 4-component current by a 
[4x4] interface conductance matrix: [Ic, Iz^ Ixi ly]^ = 
[G]4x4 [AK, AI4, AVyf. The four-components 
represent the charge 'c' and the three spin components 
X and y. Figure |2] is plotted using (a) and (b) along with 
LLG (see also supplementary section S.IB), which shows 
the output voltage per unit current and agrees well with 
experimental data^. To obtain a close match, polariza- 
tion of the magnets was adjusted to a value of 0.5, which 
is in a reasonable range^. (S.1B,C provide the parame- 
ters used). 

The spin-transport model in Fig{l]3 illustrates three 
basic ingredients: (I) the lead-nanomagnet-channel in- 
terface: G*^^, (II) nonmagnetic- lead conductance: G^^^^ 
and (III) the channel: G^^. For an interface between a 
non-magnetic channel and a Ferromagnet pointing along 

(1) can be modeled as a conductance matrix whose 
components can be written in terms of the scattering 
matrix between the plane 'C inside the channel and the 
plane 'L' inside the lead. It has been shown^ that: 
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where ^ = 2 - r^r^ - r^r* , gP = Vdr^ - r^rl , F = 1 - r^r* 
and P is the polarization. Vu^Vd being the reflection co- 
efficients for up and down spins respectively as seen from 
the plane 'C inside the channel. This is for a single con- 
duction mode. All modes have to be added together to 
represent the interface area and materials used (S.IC). 

For (II) i.e. G^^^^ we construct the full conductance 
matrix for the contact by placing the interface conduc- 
tance in series with a Il-conductance network whose se- 
ries (se) and shunt (sh) components are given by (p: re- 
sistivity, t. length , A: cross-sectional area. A: spin-flip 
length): 
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where Qse = {A/ pX) csch {£/ X) and gsh = 
(A/pA) tanh (-^/2A). These conductance matrices 
are obtained by solving standard spin diffusion equa- 
tions which are summarized in S.IA and contain all the 
physics of spin diffusion in one dimension noting that in 
non-magnetic materials there is no distinction between 
x^y and z components. The final ingredient (III) is G^^, 
which we assume to be non-magnetic in this paper and 
therefore it is adequately described using same matrices 
as G^^^^. 

The spin components Ins of the current [In] {Fig^p) 
into magnetic contact 'n' provide the spin-torque that 
enters the LLG equation describing the dynamics of 
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FIG. 2: Calculated spin- valve signal vs. input current closely 
matches the experimental results in Ref.[7] 
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where q is the charge of electron, 7 is the gyromag- 
netic ratio, a is the Gilbert damping parameter and 
Ns = Ms^/pb is the net number of Bohr magnetons 
comprising the nanomagnet {Ms = saturation magneti- 
zation and Q = volume). H = Hkz — Hdy represents 
the internal 'uniaxial anisotropy' and 'out-of-plane de- 
magnetizing' effective fields acting on the magnet. The 
last term is the spin-torque current Ist- 

Figure [3] is plotted using the self consistent model just 
described where the component of output magnetization 
along its easy axis {z) is shown. Here, we ignore the 
transit time of carriers because it is on the order of a pi- 
cosecond much less than the nanomagnet dynamics (see 
S.2). The inset shows a magnet that has ~ 5 times lower 
anisotropy field (Hk) and 7 times higher Ns (Here Q and 
not Ms has been varied to change Ns). It is clear that 
with higher Hk and lower Ns^ magnet switches faster 
but this is not in the expense of higher current since 
Ns has been reduced and it is generally accepted that 
switching current is proportional to Ns after the work of 
Sun^. Also, we find that with identical input and output 
magnets, the switching process is non-re ciproc al without 
the need to use Bennett clocking scheme^^*^^ because of 
different voltages on the two magnets (e.g. floating out- 
put in Figjl^). This causes the output magnet to switch 
faster than the input magnet. We note however that the 
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FIG. 3: For the same current I, a magnet with smaller Ns 
(7e5) is predicted to switch faster than one with Ng — 5e6. 
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FIG. 4: Various currents throughout a switching event. 

floating voltage may not be ideal for cascading circuits; 
such issues will be discussed elsewhere. Here we focus 
more on the switching energy-delay. 

With this in mind, it is instructive to look at how the 
different current components vary with time (FigQ. It 
is evident that while the charge and spin currents, Ic 
and Is continue to flow as long as a voltage is present, 
the spin torque current 1st (cf. Eqj3| that enters the 
LLG equation is time-limited: it flows only during the 
time that the magnet is switching. Indeed we find that 
Jq Igtdt = /i {2qNs). The factor /i is exactly 1 if only 
a uniaxial field is present^^ as we might have expected 
from angular momentum conservation^ (see also S.3). 
However we find that /i can be less or more than 1 when 
fields other than uniaxial are involved (S.3). From Fig|4] 
it is evident that the total charge Qtot in Eqs{l]will be 
larger than J^^"^ Igt^t and can be written as 

Qtot = Idt = ^/2/i {2qNs) (4) 

Jo 

where / is the charge current. Is is the time- average spin 
current and /2 = J Isdt/ f Igtdt is a factor reflecting the 
fact that the spin current is somewhat larger than 
\Ist\ that enters Eqjs] 

Evidently, the switching energy-delay, among other 
things, can be improved by lowering Ns. While re- 
ducing Ns, thermal stability has to be ensured which 



is determined by the activation barrier (£^5 = KuQ, 
Ku is the effective uniaxial anisotropy constant) of a 
magnet. E^j has to be ~ 10 's of kT to sustain non- 
volatility. Ns is related to through the following 
equation: Ns = Ms^/fiB = 2Ku^/ /j^bHr' Taking sta- 
bility into account, ultimate scaling requires magnetic 
materials (see e.g. Weller et al.^^) with high anisotropy 
fields {Hk = 2Ku/Ms) where only several thousand Bohr 
magnetons {Ns) can collectively give rise to stable mag- 
nets. In short, lowering Ns will keep the switching cur- 
rent low and high Hk will decrease the switching delay. 
Althoug h mak ing devices from high anisotropy magnetic 
material j^^l^^* could have experimental challenges, we be- 
lieve the underlying potential impact on lowering energy- 
delay of spin-torque switching remains valid. On the 
other hand with magnetic field switching used in schemes 
such as MQCA^, higher Hk would require a higher 
switching field; hence higher switching current. Similarly, 
the scaling of switching energy-delay based on Hk and 
Ns is not favorable to multiferroic switching of magnet 
either since higher Ku and/or lower delay would require 
higher switching voltages. Note that presently the low 
voltage operation of ASL devices is offset by the large 
total charge {Qtot ^ 2.4e7q in FigQ, arising from a com- 
bination of large magnets and low switching efficiency. If 
these numbers (i.e. Ns and /1/2^/^s) can be reduced, the 
advantages of low voltage operation (less parasitic capac- 
itance and stray charge) and non-volatility (less leakage) 
would make ASL look attractive. 

In summary, we have presented a model that combines 
the physics of spin transport with that of nanomagnet dy- 
namics that agrees well with available experimental data. 
Using this model we investigate the switching of ASL de- 
vices and show how the energy-delay scales with Ns (Eqs. 
1 and 4). It is also shown that for identical input /output 
magnets, switching can be non-reciprocal based on the 
applied voltages. Suitable cascading schemes will be dis- 
cussed elsewhere. 
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experimental data. 
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Supplementary Information 



I. MODEL 

A. Spin Diffusion Equations 

Spin diffusion equations for 1-dimensional transport 
can be written in terms of distributed resistances r^^, 
rd and spin-flip conductances gsf as: 



dj/j^u/ - q) 

dx 



d{iid/ - q) 
dx 
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where (luilJ^d are the quasi- Fermi levels for up and down 
spins. Defining the charge c and z-components of volt- 
age and current as: Vc = (/in + Md) /(-2g), h = lu^h 
and Vz ifiu - l-^d) I (-2<7), 4 = 4 - -^d,we can rewrite 



the spin diffusion equations (|5a[) and (5b) as: 
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The conductance matrices and G^^ of the main 



paper are obtained by solving Eqs.(6a) for a section 



of length i with specified values of AVc and AVz 
and then calculating Ic and 1^ at either end from 
Eqs.( 6b|6c ). The conductance matrices contain all 
the physics of spin diffusion in one dimension with 
A related to the distributed resistance and spin-flip 



conductance Qsf by A 



The advantage of 



iru+rd)gsf ' 

using charge and spin (rather than u/d) components is 
that the conductance matrices are trivially extended 
to 4-component quantities noting that in non-magnetic 
materials there is no distinction between x^y and z 
components. 



B. Analytical Calculations 

We would like to point out that the agreement (Fig. 2 
of the main paper) with experiment (Fig. 3d of Ref.[l] 
same as Ref.[7] of the main paper) could be expected 
from a simple steady state spin-transport calculation, 
as noted in the experimental papers. However, the real 
value of the full model described in the paper is twofold: 

• Our self consistent model couples spin transport 
and temporal dynamics of the magnet, thus en- 
abling us to extract the switching delay. Knowl- 
edge of delay is essential to the estimation of 
switching energy and further, to establish the scal- 
ing with Ng. However, the delay cannot be ob- 
tained from steady state spin-transport calcula- 
tions. 

• Without coupling the magnet dynamics to spin 
transport, it is impossible to demonstrate the pos- 
sibility of non-reciprocal switching as was noted in 
the paper. The analysis of any logic switch based 
on the ASL device needs to solve the spin trans- 
port model with the dynamics of input-output ma- 
gents self-consistently to analyze non-reciprocity. 

Steady state calculation: Here we outline simple esti- 
mates for the steady state calculations. In Fig 2, the 
two measured values of the non-local resistance (out- 
put voltage per unit current) correspond to collinear 
(parallel and anti-parallel) configuration of the mag- 
nets for which one can obtain an analytical expression 
using the spin-diffusion equations presented in the pre- 
vious section: 
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Pf{n) and Xf{n) are the resistivities and the spin dif- 
fusion lengths of the Ferromagnet (channel). Rsn = 
Pn^n/^n and i^i,i^2 = Pf^f/^f are the resistances 
over one spin-flip length for the channel and the input 
and output Ferromagnets respectively. PI (2) is the po- 
larization of the input (output) magnet. Equation 1 is 
in agreement with the results obtained by other authors 
(see e.g., the analysis of lateral spin valves by Taka- 
hashi and Maekawa^). Using the parameters listed in 
Table 1, we calculate the non-local resistance to change 
by ~ AmQ as observed experimentally and computed 
from our model. 

The critical current at which the output magnet 
switches can also be estimated analytically as noted 
in the experimental paper. For collinear magnets we 
can write the spin current in the output magnet as 



he (l + xi)(l + X2)-e-2^A 



(9) 
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TABLE I: List of material and geometric parameters for non-local spin valve structure of Ref.l 





Area 
nm2 


Length Resistivity 
nm (p) Q — nm 


Spin-flip length Spin resistance 
(A)nm {Rs = p\/A) VL 


Input magnet (NiFe) 


75 X 170 


20 


170 


5 


0.07 


Input lead (Au) 


75 X 170 


50 


70 


10 


0.05 


Output magnet (NiFe) 


80 X 170 


4 


170 


5 


0.06 


Output lead (Au) 


80 X 170 


50 


70 


10 


0.05 


Channel (Cu) 


170 X 65 


^ = 270 


7 


1000 


0.63 


Left lead (Cu) 


170 X 65 ^ 


1 = 10000 


7 


1000 


0.63 


Right lead (Cu) 


170 X 65 ^ 


2 = 10000 


7 


1000 


0.63 



which is approximately 10% for the parameters in the 
experiment. For the output magnet, using the expres- 
sions obtained by Sun^, we estimate a critical spin cur- 
rent of 



2q 

Is,cr = -^01 {2KuQ ■ 



KdQ)) = 0.46mA (10) 



Here we used a Gilbert damping parameter a = 0.007 
and estimated Ku HkMs/2 and Kd 27tM^ 
using Hk = 80 Oe, Ms = 780 emu/cc. Noting that 
the input current is approximately 10 times the spin 
current delivered to the output magnet, a switching 
current of 4.6 mA is comparable to that observed in 
the experiment as well as Fig.l of our paper. Also 
note that the experiment is performed at lOK so the 
zero temperature estimations from LLG are relatively 
accurate. 



II. TRANSIT TIME 

Here we would like to estimate the transit time of car- 
riers due to the diffusion in channel from input magent 
to the output magnet. As we will see shortly, this time 
scale is much faster than nanomagnet dynamics and 
because of this, the time of flight of carriers has been 
neglected in our models presented in the paper. This 
time of flight can be estimated as: 



2D 



(11) 



Where L is the distance traveled and D is the diffusion 
coefficient. We set D = '^F^mfD- '^f = 1-57 x lO^m/s is 
the Fermi velocity for coppei^ and Xmfp is the mean 
free path which can be calculated using: 



C. Parameters used for the experimental 
benchmark of Fig. 2. 

As mentioned earlier in the main paper, our model 
describes real structures fairly well as it mimics the 
switching characteristics of the experiment of Yang 
et.al.1' (Fig. 2). The material and geometric parameters 
used to benchmark our model are summarized in Table 
1. These parameters have been obtained from experi- 
mental papers^ ^ and private communications with the 
authors. Apart from the parameters listed in the table, 
we have to point out that one also needs the value of 
the so called mixing conductances^ which give the x — y 
block of the interface conductance matrix G*"^^ men- 
tioned in our paper. These values generally are calcu- 
lated using ab inito models. We use quantities provided 
by Brataas et al.^ noting that they can be estimated 
using gmix = ^q^/^/47t^ where gq is the quantum of 
conductance, kf is the wave- vector in the channel and 
A is the area of the magnet-channel interface. Note 
that and are the matrix elements that give 
rise to the field-like term and it is generally believed 
that in such metallic structures they are negligible. For 
a precise match, the polarization of the NiFe/Cu mag- 
net channel interface has been assumed to be 0.5, well 
within the experimentally measured rang 



^mfp 



amvF 



(12) 



where cr = 5.8 x l{fVt~^m~^ is the copper conductivity, 
m = 9.11 X 10~^^kg is the free electron mass, n = 8.47 x 
/m? and g = 1.6 x 10~^^C is the electron charg^. 
Using these numbers we get Amfp ^ 38 nm. One could 
also use a more conservative estimate of Amfp ^10 nm 
for less pure copper. We can now calculate the time it 
takes for the carriers to traverse a distance of L = 100 
nm in copper: 



(lOOxlO-9)^ 



ransit — 2X^fpVF ~ 2Amfp (1.5 X lO^) 

1.7 X 10"^^sec for Amfp = 38nm 
6.4 X 10"^^5ec for Amfp = lOnm 



These quantities clearly show that the dominant char- 
acteristic times are not the transit time of carriers, but 
rather they are that of nanomagnet dynamics, there- 
fore at this time they have not been included in our 
models. 
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III. CONSERVATION OF ANGULAR 
MOMENTUM 

A. Case 1: No out-of-plane demagnetizing field. 

(Hd = 0) 

We are interested in the time-integral of LLG equation 
throughout switching: 
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-m X (^m X Ins^ dt 



For the component of magnetization along its easy axis 
{z) we have: 
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In the absence of any out-of-plane demagnetizing field, 
the precessional term (m x H) has no z-component; 
hence the first term on the right hand side vanishes. 
For the second term we have: 



rtsw 


dm 


/ a 


m X — — 
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dt _ 



dt 



m^dmy- 



mydmx = 



Each integral on the right hand side of above equation 
amounts to zero because it is an integration of a total 



differential over a closed path. Noting the left hand side 
of Eq 13 for a switching event (±z to ^z) amounts to 
±2 we have: 



/* 

^0 



\Ilt\dt 



q2N, 



B. Case 2: Non-zero Out-of-plane demagnetizing 
field. {Hd -ATTMsTUyy) 

For the case where both the in-plane uniaxial field and 
the out-of-plane demagnetizing field are also present 
the first term on the r.h.s of Eqjl3]has a non- vanishing 
component: 



m X H 



dt 



mxmydt 



Note that this term can amount to net positive or neg- 
ative values. This depends on various conditions such 
as the ratio of the uniaxial anisotropy field to the de- 
magnetizing field and the amount of current overdrive. 
Based on the simulations that we have performed, we 
see the following range: 

1'"^' \FJdt = M{2Ns), / = 0.5 -1.8 

which never the less gives a value on the order of 2Ns 
for a switching event. 
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